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Mechanical  stiffness  and  dissipation  in  ultrananocrystalline  diamond  microresonators 
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6 Department  of  Mechanical  Engineering  and  Applied  Mechanics,  University  of  Pennsylvania,  Philadelphia,  Pennsylvania  19104,  USA 

(Received  27  January  2009;  published  2  June  2009) 

We  have  characterized  mechanical  properties  of  ultrananocrystalline  diamond  (UNCD)  thin  films  grown 
using  the  hot  filament  chemical  vapor  deposition  (HFCVD)  technique  at  680  °C,  significantly  lower  than  the 
conventional  growth  temperature  of  ~800  °C.  The  films  have  ~4.3%  sp2  content  in  the  near-surface  region  as 
revealed  by  near  edge  x-ray  absorption  fine  structure  spectroscopy.  The  films,  ~1  pm  thick,  exhibit  a  net 
residual  compressive  stress  of  370  ±  1  MPa  averaged  over  the  entire  150  mm  wafer.  UNCD  microcantilever 
resonator  structures  and  overhanging  ledges  were  fabricated  using  lithography,  dry  etching,  and  wet  release 
techniques.  Overhanging  ledges  of  the  films  released  from  the  substrate  exhibited  periodic  undulations  due  to 
stress  relaxation.  This  was  used  to  determine  a  biaxial  modulus  of  838  ±2  GPa.  Resonant  excitation  and 
ring-down  measurements  in  the  kHz  frequency  range  of  the  microcantilevers  were  conducted  under  ultrahigh 
vacuum  (UHV)  conditions  in  a  customized  UHV  atomic  force  microscope  system  to  determine  Young’s 
modulus  as  well  as  mechanical  dissipation  of  cantilever  structures  at  room  temperature.  Young’s  modulus  is 
found  to  be  790  ±  30  GPa.  Based  on  these  measurements,  Poisson’s  ratio  is  estimated  to  be  0.057  ±  0.038.  The 
quality  factors  ( Q )  of  these  resonators  ranged  from  5000  to  16000.  These  Q  values  are  lower  than  theoretically 
expected  from  the  intrinsic  properties  of  diamond.  The  results  indicate  that  surface  and  bulk  defects  are  the 
main  contributors  to  the  observed  dissipation  in  UNCD  resonators. 

DOI:  10. 1103/PhysRevB. 79. 245403  PACS  number(s):  62.25. -g,  62.40.+i,  62. 30.+d,  62.20.dj 


I.  INTRODUCTION 

Nanocrystalline  materials  present  a  rich  set  of  opportuni¬ 
ties  to  explore  the  physics  governing  how  defects  such  as 
grain  boundaries,  impurities,  and  dangling  bonds  affect  the 
mechanical  performance.  Diamond,  the  stiffest  material 
known,  can  be  grown  as  a  conformal  thin  film  on  a  variety  of 
surfaces  in  a  form  known  as  ultrananocrystalline  diamond 
(UNCD)  (grain  size  —2-5  nm).1"3  UNCD  films  were  origi¬ 
nally  grown  using  microwave  plasma  chemical  vapor  depo¬ 
sition  (MPCVD),  and  Young’s  modulus4,5  and  hardness5  of 
films  grown  using  this  well-established  recipe  have  been 
measured.  These  films  exhibit  smooth  surfaces1  and  me¬ 
chanical  properties  close  to  single-crystal  diamond  despite 
the  high  volume  fraction  of  grain  boundaries  and  the  pres¬ 
ence  of  atoms  such  as  hydrogen.6  However,  Poisson’s  ratio 
of  UNCD  has  never  been  reported,  and  measurements  of 
mechanical  dissipation  (quality  factor)  are  scant.  While  mod¬ 
est  deviations  in  the  elastic  properties  (modulus,  Poisson’s 
ratio)  from  the  values  for  bulk  diamond  can  be  expected  due 
to  the  polycrystalline  nature  of  UNCD  or  the  possibility  of 
defects,  the  quality  factor  is  far  more  difficult  to  predict  since 
mechanical  dissipation  has  a  large  number  of  possible 
sources  which  are  often  very  sensitive  functions  of  several 
material  parameters. 

In  general,  UNCD  and  nanocrystalline  diamond  (grain 
size  —10  nm-100  nm)  thin  films  exhibit  a  unique  combi¬ 
nation  of  properties  suitable  for  applications  to 
microelectromechanical/nanoelectromechanical  system 
(MEMS/NEMS)  devices,2  such  as  radio  frequency  MEMS/ 
NEMS  resonators. 7-11  This  is  due  to  their  superior  physical 
properties  such  as  a  high  Young’s  modulus,  their  stable  and 


inert  surface  with  low  adhesion  to  other  materials,  and  their 
good  tribological  performance. 12-14  High-frequency  devices 
based  on  MEMS  or  NEMS  structures  have  a  large  surface  to 
volume  ratio  and  can  have  high  dissipation  if  surface  effects 
dominate  the  loss  mechanisms.1516  However,  the  surface  sta¬ 
bility  and  high  stiffness  of  UNCD  films  suggest  it  has  poten¬ 
tial  for  high-frequency  resonators  with  high  quality  factors. 

The  ability  to  grow  uniform  and  smooth  UNCD  films 
over  a  large  area  at  relatively  low  temperatures  (<500  °C) 
has  enabled  the  cointegration  of  complementary  metal- 
oxide-semiconductor  (CMOS)  electronics  with  UNCD 
MEMS  devices  (e.g.,  resonators  and  switches).1,2  The  me¬ 
chanical  properties  of  UNCD  films  grown  at  low  tempera¬ 
tures  play  an  important  role  in  the  suitability  of  these  films 
for  applications  involving  monolithic  integration  of  MEMS/ 
NEMS  and  CMOS  devices.  In  this  paper,  we  present  studies 
of  Young’s  modulus,  Poisson’s  ratio,  and  mechanical  dissi¬ 
pation  of  UNCD  films  grown  at  680  °C  using  the  hot  fila¬ 
ment  chemical  vapor  deposition  (HFCVD)  technique.  We  de¬ 
termined  the  Young’s  modulus  of  these  films  by  fabricating 
UNCD  microcantilevers  and  measuring  the  resonance  fre¬ 
quency  of  the  fundamental  flexural  mode  of  the  cantilevers. 
Underetching  of  the  UNCD  films  to  produce  overhanging 
ledges  resulted  in  periodic  wrinkles  due  to  residual  compres¬ 
sive  stress  in  the  film.  The  amplitude  and  wavelength  of  the 
fully  relaxed  overhanging  portion  were  measured  to  deter¬ 
mine  the  biaxial  modulus  of  the  film  and,  by  comparison 
with  the  Young’s  modulus,  Poisson’s  ratio.  Ring-down  mea¬ 
surements  of  resonantly  excited  cantilevers  were  used  to  de¬ 
termine  mechanical  dissipation  in  the  UNCD  at  —kHz  fre¬ 
quencies. 
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II.  DISSIPATION  IN  MICROCANTILEVERS 


Dissipation  in  resonators  can  be  intrinsic  or  extrinsic  in 
origin.  Viscous  damping  from  a  surrounding  fluid  medium 
(e.g.,  air),  clamping  losses,  or  dissipation  due  to  the  measur¬ 
ing  system  itself  are  examples  of  extrinsic  losses.  Under 
UHV  conditions,  air  damping  is  eliminated.17  A  lower  bound 
on  dissipation  from  the  measuring  system  can  be  determined 
by  carrying  out  a  measurement  on  a  high  Q  structure. 

There  are  several  models18^20  for  quality  factors  limited 
by  clamping  losses,  represented  as  GciamPing>  for  out  of  plane 
flexural  cantilever  resonators  attached  to  supports  of  differ¬ 
ent  dimensions.  Jimbo  et  al.20  predicted  the  quality  factor 
limited  by  clamping  loss  of  cantilever  beam  of  infinite  width 
attached  to  semi-infinite  base  to  be  given  by 

Gclamping^pj  >  (1) 

where  L  is  the  length  and  t  is  the  thickness  of  the  cantilever. 
Later  Cross  et  al. 17  proposed  equations  for  out  of  plane  flex¬ 
ural  resonators  of  finite  width  attached  to  thickness-matching 
bases.  For  the  base  thickness  ( tb )  smaller  than  the  wave¬ 
length  of  the  elastic  wave  in  the  base  (\b),  quality  factors 
limited  by  clamping  losses  to  the  thickness -matching  base 
are  proportional  to  the  aspect  ratio  of  the  levers  (L/w). 18,19 
However  this  model  does  not  take  into  account  the  fact  that 
thickness-matching  bases  can  be  overhangs  where  the  energy 
can  be  reflected  back  to  the  oscillator.18  Photiadis  et  al ,19 
developed  a  model  for  clamping  losses  of  flexural  cantilever 
resonators  with  finite  (tb<\b,tb>t)  and  semi-infinite  (th 
>X.J  bases.  For  finite  ( tb<\b )  base,19 

Gclamping  ^  •  (2) 


Intrinsic  dissipation  processes  relevant  to  UNCD  cantile¬ 
vers  include  thermoelastic  dissipation  (TED),2 1-2 1  phonon- 
phonon  dissipation,24  and  strain-assisted  relaxation  of 
defects.25  Thermoelastic  dissipation  occurs  mainly  due  to  the 
temperature  gradients  set  up  by  the  oscillations  in  the  solid 
which  causes  local  volume  changes.  This  results  in  irrevers¬ 
ible  heat  flow  across  the  temperature  gradient  resulting  in 
dissipation  of  the  mechanical  energy.  Zener  et  al .2I  devel¬ 
oped  an  approximate  equation  for  thermoelastic  dissipation, 
which  for  flexural  beams  is  given  by 


G"1 


EorT  <otz 
pCP  1  +  (wr.)2’ 


(3) 


where  a,  CP,  and  T  are  the  coefficient  of  thermal  expansion, 
specific  heat,  and  temperature  of  the  material,  respectively,  a> 
is  the  oscillation  frequency,  and  r:  is  the  relaxation  time  for 
TED  given  by 


t 


2 


Tz  ~  ttD  ’ 


(4) 


where  D  is  the  thermal  diffusivity.  For  UNCD,  D  will  be 
much  smaller  than  single-crystal  diamond  due  to  the  lower 
thermal  conductivity  of  UNCD26  caused  by  the  large  fraction 
of  grain  boundaries  in  UNCD.  Precise  solutions  for  the  ther¬ 


moelastic  dissipation  in  thin  beams  were  developed  by  Lif- 
shitz  et  al  22  However,  TED  calculated  from  an  exact  solu¬ 
tion  derived  by  Lifshitz  and  an  approximate  solution  given 
by  Zener  are  in  close  agreement.  TED  in  polycrystalline  ma¬ 
terials  is  discussed  by  Srikar  et  al 23  For  low-frequency 
UNCD  microcantilevers,  TED  due  to  Zener  damping  will  be 
much  larger  than  TED  due  to  intracrystalline  and  intercrys¬ 
talline  damping  mechanisms  which  are  active  at  GHz  fre¬ 
quencies.  This  is  due  to  the  extremely  high  thermal  conduc¬ 
tivity  of  individual  diamond  grains  and  their  nanoscale 
dimension.2' 

Another  important  form  of  intrinsic  dissipation  occurs  due 
to  the  interaction  between  thermal  phonons  and  the  mechani¬ 
cal  oscillations  of  the  resonator.  If  the  wavelength  of  the 
acoustic  wave  is  much  larger  than  the  mean  free  path  of  the 
thermal  phonons,  then  this  low-frequency  mechanical  oscil¬ 
lation  locally  perturbs  the  thermal  phonon  distributions  away 
from  equilibrium.25  Restoring  this  equilibrium  cause  dissipa¬ 
tion  in  the  oscillating  beam.  Dissipation  by  such  an  interac¬ 
tion  is  given  by 

_  CTr  (orph 

'ph  pv2  1  +  («Tph)2’ 

where  yis  Griineisen’s  constant,  v  is  the  sound  velocity,  C  is 
the  heat  capacity  per  unit  volume,  and  rph  is  the  phonon 
relaxation  time,  which  for  flexural  beams  is  given  by 


Gph 


Tph: 


3/c 

Cv2n 


(6) 


where  vD  is  the  Debye  sound  velocity,  which  depends  of  the 
transverse  and  longitudinal  sound  velocities,25  and  k  is  the 
thermal  conductivity. 

However,  dissipation  in  most  oscillators,  including  canti¬ 
levers,  is  usually  dominated  by  intrinsic  dissipation  due  to 
the  stress-assisted  relaxation  of  defects.  When  subjected  to 
stress,  impurity  atoms,  dislocations,  and  other  defects  at 
grain  boundaries  and  surfaces  undergo  a  transition  from  one 
state  to  another,  similar  to  two  level  systems.  Dissipation  in 
such  systems  exhibit  Debye  peaks  with  a  characteristic  re¬ 
laxation  time  t  which  is  unique  to  that  transition.  The 
frequency-dependent  quality  factor  limited  by  such  transi¬ 
tions  is  given  by27 


G  Defect  '4 


1  +  (<ur)2 


(7) 


where  A  is  a  constant  which  depends  on  the  nature  of  the 
defect  and  defect  concentration  (number  of  defects  per  unit 
volume)  and  r  is  the  relaxation  time  which  follows  the 
Arrhenius  relation  given  by27 


1 

T 


(8) 


where  EA  is  the  activation  energy  for  the  process,  kB  is  the 
Boltzmann  constant,  and  t0  is  the  characteristic  atomic  vi¬ 
bration  frequency. 

Another  important  form  of  intrinsic  dissipation  occurs  at 
interfaces  of  high-frequency  (rf)  resonators,  especially  at  the 
surfaces  of  resonators.  Dissipation  at  surfaces  could  domi- 
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FIG.  1.  NEXAFS  spectra  from  a  H-terminated  UNCD  film  and 
single-crystal  diamond  (offset  for  clarity).  UNCD  exhibits  a  sp2 
peak  at  285  eV,  a  diminished  exciton  peak  at  289  eV,  and  shallower 
band  gap  at  302  eV,  all  due  to  the  fractional  presence  of  non¬ 
diamond  bonding. 

nate  the  losses  in  oscillators  that  have  a  high  surface-to- 
volume  ratio1516  and  could  arise  from  the  presence  of  a  large 
number  of  defects  at  the  surface,  the  presence  of  extraneous 
material  (contamination),  defects  induced  by  processing 
steps  such  as  etching,  or  simply  roughness  on  the  surfaces  of 
the  structure. 


HI.  FILM  GROWTH  AND  CHARACTERIZATION 

UNCD  films  (known  as  Aqua  25  by  Advanced  Diamond 
Technologies  Inc.)  were  grown  using  the  HFCVD  technique 
at  680  °C  on  silicon  wafers  (diameter=  150  mm)  using  a 
predominantly  methane/hydrogen  growth  chemistry.  Silicon 
wafer  were  ultrasonic  ally  seeded  with  ultradispersive  detona¬ 
tion  diamond  (UDD)  powder  solution  before  the  growth.  The 
average  film  thickness  was  measured  by  the  NANOSPEC 
technique  to  be  1.002  pm,  uniform  to  within  11%  across  the 
wafer  with  the  thickness  greatest  at  the  center.  These  films 
exhibited  a  residual  compressive  stress  of  370  ±  1  MPa  mea¬ 
sured  using  a  Tencor  Flexus  2320A  stress  measurement  tool. 
NEXAFS  spectroscopy  was  used  to  determine  the  chemical 
bonding  nature  of  the  films,  especially  for  the  presence  of 
sp3  and  sp 2  content  in  the  near-surface  region  (top  — 4  nm) 
of  the  films.25  A  NEXAFS  spectrum  taken  on  the  topside  of 
a  hydrogen-terminated  UNCD  film  grown  with  the  condi¬ 
tions  indicated  above  is  shown  in  Fig.  1  (black  curve),  where 
a  NEXAFS  spectmm  from  a  single-crystal  diamond  is  shown 
for  comparison  (gray  curve).  Characteristic  diamond  features 
observed  in  Fig.  1  include  the  exciton  peak  at  289.3  eV  and 
second  band  gap  at  —302  eV  characteristic  of  diamond.  A 
small  amount  of  non-diamond  bonding  in  the  UNCD  film  is 
revealed  by  the  presence  of  ,v/;2-honded  carbon  correspond¬ 
ing  to  the  peak  in  the  electron  yield  centered  at  285.0  eV. 
Based  on  NEXAFS  measurements  on  HFCVD  grown  UNCD 
films  and  the  corresponding  reference  spectra  obtained  on 
single-crystal  diamond  and  highly  oriented  pyrolytic  graphite 
(HOPG),  we  estimate  the  sp2  content  of  the  films  to  be  5.8% 
using  the  normalization  method  described  by  Lenardi  el  a l.23 
The  sp2  content  observed  in  the  HFCVD  UNCD  films  stud¬ 


FIG.  2.  (Color  online)  (a)  SEM  image  of  unreleased  UNCD 
cantilever  die  which  features  four  different  cantilevers  with  lengths 
ranging  from  100  to  400  pm,  with  a  thickness  of  1  pm.  (b)  Opti¬ 
cal  microscope  image  of  the  back-cleaved  cantilever  die. 

ied  here  is  similar  to  the  2-5  %  sp2  content  observed  in 
MPCVD  films  grown  at  800  °C.3  Most  of  this  ,v/?2-bondcd 
carbon  is  present  at  the  grain  boundaries;  some  fraction  may 
be  due  to  adsorbates  and  surface  reconstruction.  Hydrogen 
termination29  of  the  surface  using  a  hot  filament  process  re¬ 
sulted  in  the  reduction  in  sp2  content  in  the  near-surface 
region  to  4.3%.  Thus,  this  percentage  represents  a  likely 
lower  bound  for  the  sp2  content  of  the  bulk  of  the  film  as¬ 
suming  a  uniform  film  structure  with  thickness.  The  UNCD 
films  studied  here  exhibit  an  RMS  roughness  of  —10  nm 
over  10  X  10  pm2  area  measured  by  tapping  mode  AFM  im¬ 
aging. 

IV.  FABRICATION  AND  EXPERIMENTAL  SETUP 

For  cantilever  fabrication,  a  thick  sputtered  oxide  was  de¬ 
posited  on  the  UNCD  surface,  and  subsequently  patterned 
and  etched  to  serve  as  a  hard  mask  for  etching  the  UNCD 
layer  to  produce  the  cantilevers.  The  UNCD  layer  was  etched 
using  oxygen-based  plasma  in  a  reactive  ion  etching  (RIE) 
process  to  define  the  UNCD  cantilevers.  The  wafers  were 
diced  to  commercial  AFM  chip  specifications  (4.4 
X  1.6  mm2  dies)  to  enable  insertion  into  commercial  AFMs 
for  measuring  the  cantilever  properties.  Each  die  contained 
four  cantilevers  as  shown  in  the  SEM  image  in  Fig.  2(a), 
with  the  cantilever  lengths  of  100,  250,  300,  and  400  pm. 
The  length-to-width  ratio  is  7.5  for  all  the  cantilevers.  The 
oxide  hard  mask  and  sacrificial  oxide  underneath  the  canti¬ 
levers  were  subsequently  wet  etched  simultaneously,  using  a 
49%  hydrofluoric  acid  (HF)  by  volume  solution.  The  sub¬ 
strates  were  cleaved  so  that  the  cantilevers  did  not  have  any 
substrate  material  underneath  them  [Fig.  2(b)].  This  was 
done  to  ensure  that  reflected  light  was  sampled  only  from  the 
cantilever  and  not  from  the  more  reflective  Si  substrate.  Dies 
are  further  etched  with  either  a  SF6  or  a  XeF2  dry  etch  chem¬ 
istry  to  study  the  effect  of  overhang  on  mechanical  proper¬ 
ties.  The  thickness  of  each  cantilever  was  measured  by  SEM 
and  ranged  from  820  nm  to  1.01  pm.  The  minor  differences 
in  thickness  are  from  thickness  variation  of  the  film  as  de¬ 
posited  or  due  to  subsequent  processing.  Dies  were  then 
cleaned  using  a  piranha  solution  (3:1  concentrated  sulfuric 
acid  to  30%  hydrogen  peroxide  solution)  before  being  tested 
under  UHV  (2  X  10_1°  Torr)  conditions  on  a  custom-built 
AFM  stage  inserted  into  a  RHK  UHV  350  AFM.  The  AFM 
stage  consisted  of  a  piezoceramic  glued  onto  an  invar  holder 
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FIG.  3.  (Color  online)  Custom-made  AFM  stage  with  UNCD 
cantilever  die  mounted  on  a  tilted  piezoceramic.  This  stage  is  in¬ 
serted  into  a  RHK  AFM  sample  holder.  A  thermocouple  is  attached 
to  the  die  to  measure  the  temperature. 

tilted  at  22.5°  (so  as  to  reflect  incident  laser  light  off  of  the 
cantilever  to  the  photo  detector)  as  shown  in  Fig.  3  (see 
schematic).  Cantilever  dies  were  mounted  on  top  of  the  pi¬ 
ezoceramic  and  the  cantilever  deflection  was  measured  using 
the  AFM  photo  detector.  Resonant  frequencies  were  identi¬ 
fied  by  sweeping  the  input  frequency  to  the  piezo  across  the 
expected  frequency  range  for  the  cantilever  beams.  To  deter¬ 
mine  the  quality  factor,  a  ring-down  measurement  approach 
was  employed  whereby  the  excitation  was  stopped  and  the 
subsequent  cantilever  motion  recorded  using  a  Tektronix 
3014B  oscilloscope.  An  exponential  fit  to  the  ring-down 
curve  gives  a  measure  of  dissipation  in  the  fixed-free  beams. 
The  sample  location  was  used  instead  of  the  regular  cantile¬ 
ver  location  as  it  allowed  us  to  independently  cool  the  can¬ 
tilever  using  the  variable  temperature  sample  stage  of  the 
AFM.  Results  at  varying  temperatures  will  be  reported  sepa¬ 
rately. 

Dies  without  cantilevers  were  also  studied  to  determine 
the  biaxial  modulus.  XeF2  dry  etching  chemistry  is  used  to 
produce  an  uninterrupted  overhang  of  UNCD  films  grown  on 
silicon.  The  amplitude  and  wavelength  of  the  overhanging 
portion  were  measured  using  a  Zygo  white  light  scanning 
interferometer  (New View  6K).  Different  overhangs  were 
produced  until  the  film  is  completely  relaxed  to  determine 
the  biaxial  modulus. 

V.  RESULTS  and  DISCUSSION 
A.  Young’s  modulus 

The  measured  resonant  frequencies  of  cantilevers  of  dif¬ 
ferent  lengths  and  overhangs  are  plotted  in  Fig.  4.  The  reso¬ 
nant  frequency  for  the  nth  mode  of  an  undamped  freely  vi¬ 
brating  cantilever  is  given  by27 


where  L  and  t  are  the  length  and  thickness  of  the  beam, 
respectively,  and  E  and  p  are  Young’s  modulus  and  density, 
respectively.  (in  is  a  constant  equal  to  —1.875,  4.694,  and 
7.855  for  n=  1,2,3,  respectively.  Therefore,  the  Young’s 
modulus  of  the  cantilever  material  can  be  determined  if  the 
mass  density  and  the  dimensions  are  known.  Measurements 
of  resonance  frequencies  vs  1  /  L1  for  all  cantilevers  with 
overhangs  ranging  from  5  to  85  /rm  are  shown  in  Fig.  4(a). 


FIG.  4.  (Color  online)  (a)  First  flexural  resonance  frequency  of 
the  cantilevers  as  a  function  of  the  inverse  square  of  the  cantilever 
length,  (b)  A  comparison  of  the  modulus  determined  from  Eq.  (9)  as 
a  function  of  actual  cantilever  length,  before  and  after  adding  an 
effective  length  of  36  pm  to  the  actual  length  of  the  levers.  All 
three  cantilevers  were  on  a  single  die  with  an  85  pm  overhang. 

In  principle,  if  the  mass  density  is  specified,  Eq.  (9)  could 
then  be  used  to  determine  Young’s  modulus.  However,  the 
presence  of  an  overhang  at  the  cantilever  base  alters  the  reso¬ 
nant  frequency  from  the  value  predicted  by  Eq.  (9)  by  an 
amount  which  increases  as  the  lever  gets  shorter.  Therefore, 
we  employed  two  methods  to  correct  for  this  effect.  First,  we 
determined  the  modulus  numerically  using  the  COMSOL  mul¬ 
tiphysics  package  (COMSOL  Inc.  Burlington,  MA,  USA), 
accounting  for  the  measured  overhang.  This  yielded  a 
Young’s  modulus  of  790  ±  30  GPa.  Second,  the  Young’s 
modulus  was  independently  estimated  by  adding  an  effective 
length  chosen  to  give  a  constant  Young’s  modulus  for  all  the 
levers  of  different  lengths.30  Figure  4(b)  shows  such  an  esti¬ 
mate  of  the  Young’s  modulus  for  cantilevers  with  85  pm 
overhang  obtained  by  adding  36  pm  to  the  actual  length  of 
the  cantilevers.  Based  on  this  estimate,  we  determine  the 
modulus  to  be  792  ±  38  GPa. 

Both  of  the  calculations  of  the  Young’s  modulus  described 
above  depend  linearly  on  the  density  of  the  UNCD  films  as 
shown  in  Eq.  (9).  The  density  is  not  known  with  certainty. 
For  the  calculations  presented  here,  the  density  of  the  UNCD 
is  assumed  to  be  equal  to  that  of  single-crystal  diamond 
(3500  kg/m3).  The  reported  uncertainty  of  the  Young’s 
modulus  of  individual  cantilevers  was  determined  from  the 
accuracy  of  the  thickness  of  the  beams  measured  in  SEM, 
which  leads  to  an  error  of  —4%  in  the  measured  Young’s 
modulus.  The  Hershey-Kroner-Eshelby  averaging  procedure 
for  polycrystalline  diamond  predicts  an  isotropic  Young’s 
modulus  of  1143  GPa.31  Thus,  the  measured  Young’s  modu¬ 
lus  is  —25-30%  lower  than  the  theoretically  calculated 
value  for  microcrystalline  diamond  with  randomly  oriented 
grains.31  Our  measured  value  is  —5-15%  lower  than  the 
experimentally  determined  value  for  the  modulus  of  UNCD 
films  grown  using  MPCVD  at  —800  °C.  Measurements  on 
these  traditional  UNCD  were  performed  using  nanoindenta¬ 
tion  (864  GPa)5  and  cantilever  beam  deflection  (916-959 
GPa).4  However,  MPCVD-grown  NCD  fixed-fixed  resonator 
structures  yielded  Young’s  modulus  values  ranging  from 
680-980  GPa,9-11  and  MPCVD-grown  UNCD  fixed-fixed 
resonators  grown  at  550  °C  yielded  a  modulus  of  710  GPa.2 

The  reasons  for  the  somewhat  lower  Young’s  modulus  of 
the  UNCD  films  studied  here  compared  to  MPCVD  grown 
traditional  UNCD  films  discussed  are  not  understood  at  this 
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time.  Theoretical  studies32-34  of  mechanical  properties  of 
UNCD  have  indicated  that  the  observed  Young’s  modulus 
strongly  depends  on  the  sp2  content  as  well  as  volume  frac¬ 
tion  of  grain  boundaries.  Philip  et  al 25  showed  that  NCD 
films  (columnar  grains,  grain  size  10-100  nm)  can  have 
Young’s  moduli  in  the  range  of  500-1120  GPa  depending  on 
the  initial  nucleation  density.  Thinner  films  of  UNCD  are 
reported  to  have  a  reduced  Young’s  modulus  compared  to 
single-crystal  diamond,36  possibly  due  to  the  enhanced  influ¬ 
ence  of  the  inhomogeneous  seed  layer  in  the  films  which  can 
be  geometrically  irregular  and  can  possess  enhanced  sp2 
bonding.1314  To  understand  the  origins  of  these  differences 
fully,  it  would  be  desirable  to  perform  a  systematic  series  of 
studies  comparing  UNCD  films  grown  using  HFCVD  and 
MPCVD  techniques  at  the  same  temperature,  with  similar 
chemistries,  such  that  both  films  show  similar  content  of  sp2 
and  sp 3  bonds,  and  comparable  nanostructures,  as  well  as  to 
systematically  vary  the  thickness.  However,  it  can  be  specu¬ 
lated  that  the  different  Young’s  modulus  observed  for  our 
films  compared  to  traditional  MPCVD  UNCD  films  and 
MCD  films  may  be  attributed  to  the  volume  fraction  of  grain 
boundaries  which  accommodate  some  amount  of  disordered 
carbon,  dangling  bonds,  and  hydrogen,  and  possibly  to  the 
seeding  method  used  and  the  resulting  initial  nucleation  den¬ 
sity. 


B.  Poisson’s  ratio 


The  biaxial  modulus  ( Eb )  of  the  thin  films  can  be  esti¬ 
mated  by  measuring  the  period  and  the  amplitude  of  the 
sinusoidal  curvature  of  the  fully  relaxed  overhang  of  the  film 
if  the  residual  compressive  stress  is  known.37,38  To  ensure  the 
film  is  completely  relaxed,  the  film  has  to  be  underetched 
until  the  ratio  of  amplitude  to  wavelength  does  not  change 
with  further  underetching.  For  the  fully  relaxed  thin  films 
with  sinusoidal  amplitude  A0  much  smaller  than  the  wave¬ 
length,  the  biaxial  modulus  is  given  by37,39 


E 


1  -  v 


(10) 


where  cr0  is  the  residual  compressive  stress  and  v  is  the  Pois¬ 
son’s  ratio  of  the  film.  To  eliminate  the  effects  of  boundaries, 
the  overhanging  ledge  should  be  several  wavelength  long  (in 
the  direction  parallel  to  the  edge  of  the  substrate),  and  the 
underetch  should  be  more  than  half  of  the  wavelength  (X/2) 
of  relaxed  overhang.  If  the  Young’s  modulus  and  the  residual 
compressive  stress  in  the  film  are  measured  independently, 
then  we  can  determine  the  Poisson’s  ratio  v  using 
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Residual  stress  in  the  film  was  relieved  by  undercutting 
the  substrate  using  dry  etching  as  described  above.  Figure  5 
shows  the  ratio  of  the  amplitude  of  the  undulation  A0  to  the 
period  X  of  the  overhanging  portion  of  the  film  as  a  function 
of  the  overhang  length.  For  overhangs  less  than  50  pm,  pe¬ 
riodic  undulations  are  barely  observed.  Beyond  an  overhang 
of  160  pm  (~X/2),  the  ratio  stabilizes  to  0.0067  and  thus 
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FIG.  5.  (Color  online)  Ratio  of  amplitude  to  wavelength  of  the 
overhang  as  a  function  of  the  overhang  length.  For  overhangs  be¬ 
tween  50  and  160  pm,  the  UNCD  film  is  not  fully  relaxed  from 
residual  stress  present  in  the  substrate.  Beyond  160  ^m,  A0/X  does 
not  change  with  further  etch  depths. 

the  film  is  fully  relaxed  and  the  remaining  substrate  does  not 
change  the  stress  relaxation  behavior  of  the  edge  of  the  film. 
Figure  6(a)  shows  the  profile  of  one  such  overhang  taken 
using  scanning  white  light  interferometer.  Figure  6(b)  shows 
the  line  profile  of  the  overhang  at  the  ridge  of  the  sinusoidal 
undulation.  Based  on  these  measurements,  we  estimate  the 
biaxial  modulus  to  be  838  ±  2  GPa.  Using  the  Young’s 
modulus  of  790  ±30  GPa  determined  using  the  resonance 
measurements  described  earlier,  the  calculated  Poisson’s  ra¬ 
tio  of  UNCD  is  then  0.057  ±  0.038.  The  large  relative  error  is 
due  to  the  subtraction  of  two  large,  similar  values  E  and  Eh 
to  produce  a  much  smaller  value. 

The  theoretical  value  of  Poisson’s  ratio  of  polycrystalline 
diamond  with  randomly  oriented  grains  and  ignoring  the 
contribution  of  grain  boundaries  is  0.069,  based  on  the 
Hershey-Kroner-Eshelby  averaging  method.31  A  low  Pois¬ 
son’s  ratio  is  a  unique  property  of  tetrahedral  sp 3  bonded 
carbon,  which  shows  a  higher  elastic  resistance  to  bond 
bending  compared  to  bond  stretching.  However  for  UNCD 
films,  which  have  a  significant  proportion  of  grain  bound¬ 
aries  containing  sp2  bonded  carbon,  the  Poisson’s  ratio  may 
be  expected  to  be  higher  than  that  for  theoretically  predicted 
value  for  polycrystalline  diamond.  Indeed,  the  measured 
Poisson’s  ratio  of  tetrahedral  amorphous  carbon  (ta-C)  films 
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FIG.  6.  (Color  online)  (a)  Interferometric  image  of  the  UNCD 
film  with  a  200  pm  etched  overhang,  (b)  Line  scan  of  the  section 
shown  in  7a  showing  a  near  saturation  beyond  160  pm.  The  appar¬ 
ent  jump  in  the  height  at  the  base  of  the  overhang  is  simply  due  to 
the  change  in  reflectivity  from  the  more  reflective  underlying 
substrate. 
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Dissipation  in  UNCD  cantilevers 
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FIG.  7.  (Color  online)  (a)  Quality  factor  of  UNCD  cantilevers 
with  different  cantilever  and  overhang  lengths,  (b)  Typical  ring- 
down  measurement  with  curve  fitting  equation. 


containing  much  higher  sp 2  content  is  0.202  ±  0.054.40  In 
comparison  with  polycrystalline  diamond,  UNCD  grain 
boundaries  are  atomically  abrupt.2  However,  because  of  the 
small  grain  size  of  UNCD,  a  significant  fraction  of  atoms  are 
in  the  grain  boundaries  or  very  close  to  the  grain  boundaries 
in  UNCD.  This  significantly  influences  the  bonding 
nature41-42  and  hence  the  mechanical  behavior  of  the  UNCD 
films.  The  influence  on  Poisson’s  ratio  for  this  special  case 
has  just  recently  been  modeled.  A  theoretical  predictions  for 
Poisson’s  ratio  of  UNCD  (calculated  from  Young’s  modulus 
and  the  bulk  modulus)  yielded  a  value  of  0.07  for  a  grain  size 
of  4.41  nm  with  significant  (9.12%)  non-diamond  (sp1) 
content,34  substantially  more  than  the  sp1  content  of  4.3% 
that  we  measured  for  our  films  using  NEXAFS.  This  sug¬ 
gests  that  a  low  Poisson’s  ratio  can  still  be  preserved  even 
with  significant  sp2  content.  However,  one  limitation  of  this 
model  is  that  it  does  not  consider  effects  of  the  presence  of 
hydrogen  in  grain  boundaries.6  As  well,  note  that  the  value  of 
Poisson’s  ratio,  like  the  value  for  E,  does  depend  on  the  mass 
density,  which  could  be  lower  for  UNCD  as  compared  to  the 
value  single-crystal  diamond,  which  we  have  used  in  our 
calculations.  A  lower  density  reduces  the  calculated  modulus 
and  hence  increases  the  measured  Poisson’s  ratio. 


C.  Dissipation  in  diamond 

Ring-down  measurements  were  conducted  at  the  funda¬ 
mental  flexural  resonant  frequencies  of  the  cantilevers  under 
UHV  conditions.  The  curve  fit  to  the  damped  oscillation 
yielded  quality  factors  in  the  range  of  5000-16000  for  all  the 
cantilevers  of  different  dimensions  and  overhangs,  with  no 
specific  dependence  on  the  dimension  of  the  cantilever  and 
the  overhang  length,  as  shown  in  Fig.  7(a).  Quality  factors 
were  observed  to  be  stable  with  time.  Figure  7(b)  shows  a 
typical  ring-down  measurement  with  a  decaying  exponential 
curve  fit  used  for  determining  the  dissipation  in  the  resona¬ 
tors.  We  have  measured  quality  factors  up  to  170000  on 
single-crystal  silicon  cantilevers  of  similar  dimension  on  the 
same  platform,  which  indicates  that  the  dissipation  due  to  the 
measuring  technique  is  minimal.  The  absence  of  any  metal 
coating  on  these  cantilevers  means  that  the  observed  dissipa¬ 
tion  can  be  attributed  to  clamping  losses  and/or  intrinsic  dis¬ 
sipation  in  the  UNCD  itself.  With  reference  to  Sec.  II,  clamp¬ 
ing  losses  for  our  cantilever  beams  with  overhangs  are 
difficult  to  model.  The  low  resonant  frequencies  (kHz)  of  our 


FIG.  8.  Solid  lines:  calculated  dissipation  in  UNCD  cantilevers 
(thickness  =  1  pm)  due  to  TED,  phonon-phonon  dissipation  as  a 
function  of  frequency.  Data  points:  measured  dissipation  at  funda¬ 
mental  flexural  resonance  of  the  cantilevers  used  in  this  study.  Since 
diamond  remains  an  ideal  Griineisen’s  solid  independent  of  the 
pressure,  a  Griineisen’s  constant  of  1  has  been  used  in  the  models. 

beams  and  high  acoustic  velocity  of  diamond  result  in  the 
wavelength  of  the  elastic  wave  in  the  base  \h  being  much 
higher  than  the  thickness  of  the  base  tb.  We  have  thickness¬ 
matching  overhangs  at  the  cantilever  base  (t=tb)  and  this 
thin  overhanging  membrane  attaches  to  thick  base  (UNCD 
deposited  on  top  of  Si).  A  more  elaborate  approach  which 
considers  the  finite  extent  of  our  overhang  and  presence  of 
different  layers  attached  to  this  thin  overhanging  membrane 
are  required  to  estimate  the  clamping  losses.  As  mentioned 
earlier,  we  did  not  see  a  systematic  dependence  of  quality 
factor  on  either  the  dimension  of  the  beam  or  the  length  of 
the  overhang.  This  strongly  indicates  that  the  losses  in  our 
beams  are  dominated  by  intrinsic  dissipation  mechanisms 
such  as  TED,  internal  friction  due  to  defects,  or  dissipation 
due  to  interaction  between  thermal  phonon  and  acoustic 
phonons. 

The  observed  dissipation  is  at  least  1000  times  higher 
than  the  theoretical  prediction  for  TED,  phonon-phonon  dis¬ 
sipation  put  together.  The  quality  factor  Q  vs  frequency  for 
UNCD  cantilevers,  jointly  with  Q  calculated  for  possible  loss 
mechanisms,  are  all  shown  in  Fig.  8.  Table  I  shows  the  con¬ 
stants  used  for  the  estimation  of  dissipation  in  UNCD  canti¬ 
levers  that  are  plotted  in  Fig.  8.  The  comparison  of  (Q's 
shown  in  Fig.  8  indicates  that  the  dissipation  is  mainly  due  to 
the  relaxation  of  large  number  of  defects  in  the  bulk  or  the 
surface  of  the  film.  Sepulveda  el  al.14  reported  a  quality  fac¬ 
tor  Q  for  polycrystalline  diamond  (PCD)  (grain  size 
—  300  nm)  fixed-free  resonators  in  the  range  of 
4000-100000,  which  includes  values  higher  than  the  Q’ s 
observed  for  UNCD  cantilevers  with  comparable  dimensions 
reported  in  this  paper.  They  suggested  that,  as  the  film  per¬ 
centage  occupied  by  the  nucleation  layer  containing  fine 
grained  diamond  is  increased,  the  quality  factor  of  the  reso¬ 
nators  is  reduced.43  This  is  consistent  with  the  idea  that  the 
higher  dissipation  in  the  UNCD  resonators  described  here  is 
mainly  due  to  the  presence  of  the  higher  proportion  of  grain 
boundaries  and  defects.  Hutchinson  et  al.10  measured  the 
temperature  dependence  of  dissipation  of  composite  metal/ 
nanocrystalline  diamond  (NCD)  fixed-fixed  beams  (colum- 
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tures  (<10  K)  indicated  the  presence  of  TLS  processes. 
Table  II  summarizes  the  results  reported  by  other 
groups9"1  f24"44’45  as  compared  to  ours  for  different  flexural 
beams  fabricated  using  NCD,  UNCD,  or  ta-C  films. 

The  temperature  dependence  of  the  dissipation  in  our 
fixed-free  beams  will  clarify  the  nature  of  dissipation  in 
UNCD  cantilever  beams.  This  will  be  discussed  in  a  separate 
publication.  However,  preliminary  investigations  point  to  the 
significant  role  of  defects  at  grain  boundaries  and  surfaces, 
indicating  that  the  variation  in  Q  that  we  observe  may  be  due 
to  variations  in  the  grain  boundary  or  surface  structure.  This 
also  suggests  that  controlling  the  grain  boundaries  and  sur¬ 
faces  may  be  the  key  to  controlling  and  increasing  the  Q 
value. 

VI.  CONCLUSIONS 

The  Young’s  modulus  of  medium-temperature  (680  °C) 
HFCVD-grown  UNCD  films  is  measured  to  be 
—790 ±30  GPa.  This  value  is  approximately  25-30  % 
lower  than  the  theoretically  predicted  value  for  the  Young’s 
modulus  of  randomly  oriented  polycrystalline  diamond  hav¬ 
ing  much  larger  grains  and  —15%  lower  than  the  experimen¬ 
tally  determined  value  for  the  Young’s  modulus  of  traditional 
UNCD  films.4  However,  it  is  significantly  higher  than  values 
measured  for  metal-coated  UNCD  beams  grown  at  lower 
temperature.2  Non-diamond  imperfections  in  the  carbon 
atom  bonding  (disordered  carbon,  dangling  bonds,  hydrogen, 
and  increased  sp2  bonding)  at  the  grain  boundaries  is  the 
likely  cause  of  the  reduction  in  the  modulus  in  UNCD  films. 
The  biaxial  Young’s  modulus  of  the  films  was  determined  to 
be  838  ±  2  GPa  by  measuring  the  amplitude  and  wavelength 
of  fully  relaxed  substrate-free  overhangs  of  the  compres- 
sively  stressed  films.  Based  on  these  measurements,  the  Pois¬ 
son’s  ratio  was  estimated  for  the  first  time  for  any  UNCD 
films  to  be  0.057  ±  0.038.  Dissipation  in  the  UNCD  cantile¬ 
vers  is  determined  using  ring-down  measurements  under 
UHV  conditions  and  ranged  from  5000-16000  at  —kHz 
resonance  frequencies.  The  UNCD  cantilever  resonators  ex¬ 
hibited  higher  dissipation  compared  to  PCD  cantilevers  with 
comparable  resonant  frequencies,  but  less  dissipation  than 


TABLE  II.  Room  temperature  quality  factor  of  flexural  resonators  from  our  work,  as  well  as  reported  by 
other  groups. 


Structure 

Resonant 

frequency 

Material 

Q 

Growth 

technique 

Actuation 

method 

Fixed-fixed“ 

17-66  MHz 

NCD 

600-2400 

MPCVD 

Magnetomotive 

Fixed-free6 

8-50  KHz 

PCD 

4000-100000 

MPCVD 

Piezoelectric 

Doubly  clamped  paddle“ 

6-30  MHz 

NCD 

2400-3500 

MPCVD 

Piezoelectric 

Fixed-free  (our  work) 

12-45  KHz 

UNCD 

5000-16000 

HFCVD 

Piezoelectric 

Fixed-fixed4 

14-157  MHz 

NCD 

-3000 

MPCVD 

Magnetomotive 

Fixed-free“ 

KHz 

ta-C 

3500 

PLD 

Piezoelectric 

“Reference  44. 
bReference  24. 
“Reference  9. 
dReference  10. 
“Reference  45. 


TABLE  I.  Constants  used  for  estimating  dissipation  due  to 
phonon-phonon  dissipation,  and  thermoelastic  dissipation. 


Parameter 

Values  used 

Cp  (J  Kg"1  K->) 

502“ 

a  (X106  K“‘) 

1“ 

y 

1“ 

T  (K) 

300 

C  (X 106J  m-3  K"1) 

1.767“ 

p  (Kg/m3) 

3500“ 

E  (GPa) 

800 

k  (W  m"1  K"1) 

10 

11.357; 

vD\  v  (X103  m  s_1) 

15.33 

“Values  used  are  for  single-crystal  diamond. 


nar  grains,  grain  size  5-15  nm).  These  high-frequency 
(MHz)  resonators  (dimensions  — 0.8  /tmXO.5  /tm 
X  10  pm)  showed  the  presence  of  a  Debye  peak  at  55  K  and 
a  dramatic  increase  in  dissipation  above  100  K.  However, 
they  reported  a  rather  lower  quality  factor  (—3000,  at  room 
temperature)  than  the  low-frequency  UNCD  cantilever  reso¬ 
nators  studied  here.  Higher  dissipation  in  the  NCD  resona¬ 
tors  can  be  attributed  to  the  fact  that  as  the  resonator  dimen¬ 
sions  become  smaller,  dissipation  due  to  surface  effects  also 
play  a  role.15,16  This  argument  is  supported  by  the  investiga¬ 
tions  of  scaling  of  dissipation  with  the  dimension  of  the 
NCD  fixed-fixed  resonators,  as  reported  by  Imboden  et  a/.44 
Factors  that  contribute  to  the  higher  dissipation  in  high- 
frequency  fixed-fixed  flexural  resonators  may  also  include 
clamping  losses,  losses  that  occur  at  metal-diamond  inter¬ 
faces  between  metal-diamond  interfaces,  and  the  dissipation 
in  the  metal  layer  of  the  composite  beams.15  Cantilever  reso¬ 
nator  structures  made  from  tetrahedral  amorphous  carbon 
(ta-C),  which  have  a  significant  sp3  content,  (80%  sp3,  20% 
sp2)  showed  a  lower  Q  of  3500.45  Dissipation  in  ta-C  indi¬ 
cated  the  processes  dominated  by  two-level  tunneling  states 
(TLS).11  Recent  investigation36  of  the  quality  factor  of  metal- 
coated  UNCD  fixed-fixed  beams  at  extremely  low  tempera¬ 
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amorphous  carbon  cantilever  resonators.45  We  attribute  dis¬ 
sipation  in  the  UNCD  resonators  mainly  to  the  presence  of 
defects  such  as  nondiamond  bonding  at  grain  boundaries  and 
surfaces. 
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